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The magnetization of quantum dots is discussed in terms of a relatively simple but exactly solv- 
able model Hamiltonian. The model predicts oscillations in spin polarization as a function of dot 
radius for a fixed electron density. These oscillations in magnetization are shown to yield distinct 
signature in the momentum density of the electron gas, suggesting the usefulness of momentum re- 
solved spectroscopies for investigating the magnetization of dot systems. We also present variational 
quantum Monte Carlo calculations on a square dot containing 12 electrons in order to gain insight 
into correlation effects on the interactions between like and unlike spins in a quantum dot. 
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I. INTRODUCTION 

As the need for nano-structures for technological ap- 
plications grows, the ability to probe and understand 
the electronic properties of these systems becomes of 
paramount importance 0, 0, Q- I n this connection, 
quantum dots (QDs), which can be viewed as artificial 
"atoms", offer unique opportunities as a nanoscale lab- 
oratory for investigating the behavior of small numbers 
of electrons and how the interplay between correlation 
and confinement effects in such systems can give rise to 
novel phenomenon such as spontaneous spin polarization 
of the electron gas 0, H, El • These and related ques- 
tions have been the subject of considerable interest in the 
recent literature 

Here we discuss how the electronic structure of QDs 
can be modeled theoretically for the purpose of gaining a 
handle on the essential phenomenology of their magnetic 
properties. The exactly solvable model Hamiltonian in- 
troduced in Ref. Q is considered first. The model of 
Ref. Q assumes a single effective interaction parame- 
ter, U, which gives the energy penalty for creating a pair 
of electrons with opposite spins. Despite its simplicity, 
this model produces considerable richness in its behavior 
and, in particular, it predicts oscillations in spin polar- 
ization with QD radius at a fixed electron density. We 
delineate the signature of spin polarization in the elec- 
tron momentum density (EMD), thus setting the stage 
for the application of momentum resolved spectroscopies 
as a window for investigating the magnetic properties of 
QDs. 

Although exactly solvable many-body models are of 
an intrinsic interest, it is important to understand the 
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nature of the parameters involved in such models via ac- 
curate first principles computations. In this connection, 
we have carried out variational quantum Monte Carlo 
(VQMC) calculations in the interacting electron gas. Re- 
sults for a 12 electron system confined within a square 
QD are presented. The computed pair correlation func- 
tions for like and unlike spins are used to deduce the ef- 
fective value of the parameter U, which enters the model 
Hamiltonian of Ref. Q . 

An outline of this article is as follows. The introduc- 
tory remarks are followed in Section II by an overview of 
the model Hamiltonian formalism of Ref. . Section III 
presents the VQMC approach and considers the example 
of a 2D square QD. Section IV presents a few concluding 
remarks. 

II. A MODEL HAMILTONIAN FOR QUANTUM 
DOTS 

Insight into properties of QDs can be obtained by con- 
sidering the relatively simple model Hamiltonian |4j 

H = Ho + Hi (1) 

where a\ a and a vcr , respectively, are the creation and an- 
nihilation operators for the one-particle state <j> va with 
eigenvalue e^. The first term (Ho) describes the noninter- 
acting system. The interaction in the second term (Hi) 
is restricted to electrons of opposite spins. The param- 
eter U here can be viewed as an average energy penalty 
for two electrons to possess opposite spins in the QD. It 
is of course energetically advantageous for electrons to 
possess the same spin because that allows the Coulomb 
energy to be lowered as the electrons are kept apart by 
the Pauli exclusion principle. 
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FIG. 1: Spin polarization per particle, £ = (iVj — N{)/N, as 
a function of the dot radius R. Solid line: interacting case; 
dashed line: weakly interacting case, U — ► 0. Simulations are 
based on the model Hamiltonian of Eq. fusing parameters 
discussed in the text. 

The model Hamiltonian of form (JIJ can be solved ex- 
actly. The solution for the many-body wavefunction has 
the form of an unrestricted Hartree-Fock wavefunction 

* = D T L> i , (2) 

where up and down arrows denote spin indices and D a = 
\<f>i >a (vj)\ is the Slater determinant formed by one particle 
states 4>i t a- Note that this solution is a spin eigenfunction 
and does not suffer from the so-called spin contamination 
problem p|. The total energy is given by 

where f ua denotes the Fermi occupation function. For a 
given total number of particles N, minimization of Eq. 
13 yields a set of nonlinear equations for the populations, 
7V| and of the up and down spins, respectively. The 
resulting splitting in energy for states of opposite spin is 

A = U(N 1 -N l ), (4) 

and it is uniform, i.e., it does not depend on the quantum 
number v. The average polarization per electron is 

C = (N r -Ni)/N. (5) 

The degeneracy between up and down spin electrons is 
thus lifted and the shell filling depends on the value of 
A, which is determined by the interaction strenght U. 

Fig. presents a few illustrative results based on the 
model Hamiltonian of Eq. ^ The details of the spe- 
cific parameters used are as follows. The non-interacting 



Hamiltonian is taken to be a 3D spherical well with po- 
tential, V(r) = -8.62 eV, for r < R, and V(r) = 0, oth- 
erwise. The electron density n in the QD is kept fixed 
corresponding to r s =5, where r s is the standard param- 
eter given by, n(47rr;?/3) = 1. The QD radius R and the 
number N of electrons are thus related via N — (R/r s ) 3 . 
r s = 5 gives a relatively low density, enabling considera- 
tion of a wide range of QD radii. The choice of U is more 
tricky since correlations in QDs are not well understood. 
However, on the basis of arguments involving a screened 
Coulomb interaction, Ref. |4j estimates U = 27.13 meV 
when N = 96 or R = 12.11 A [T^j- 

Fig-fflshows the average polarization £ as a function of 
the QD radius R for r s —5. As R increases and electrons 
are added, spin polarization ( reaches a peak each time a 
shell is half filled with up-spin electrons and falls to zero 
when the shell is completed with down-spin electrons, 
yielding a sequence of "magic numbers" , i.e., N values for 
which the QD magnetization vanishes. The oscillations in 
spin polarization are damped with increasing QD size and 
in the high R limit a paramagnetic homogeneous electron 
gas is recovered. Fig. ^ also shows that the interaction 
parameter U can give large deviations in £ from a simple 
Hund's rule filling. This is because the magnetic energy 
splitting changes with each added electron in order to 
minimize the total energy given by Eq. 

In connection with spontaneous magnetization, it is 
useful to consider the Stoner model, which is usually in- 
voked for predicting ferromagnetism in metals fll| . al- 
though it has also been applied more recently to discuss 
magnetism of nanosystems . In the Stoner model, fer- 
romagnetism results if 

ID{n) > 1, (6) 

where D(jjl) is the density of states (DOS) per unit cell 
of the spin compensated system at the Fermi level [i and 
/ is the Stoner parameter, which gives the gain in poten- 
tial energy associated with the occurrence of the ferro- 
magnetic state. Within the Density Functional Theory 
(DFT), / can be computed using the wavefunctions of 
the system at \i ^ n ^ ue case °f tuc homogeneous 

electron gas I reduces to 0] 

j = ^xc( r s) - £xc(r S )} (7) 

9(2 4 /3 _ 2) 1 > 

where e^ c and e^ c are the exchange-correlation energy 
per electron in the ferromagnetic and the paramagnetic 
electron gas, respectively. Interestingly, for the model 
Hamiltonian of Eq. ^ the connection between the aver- 
age energy penalty U for having a pair of electrons with 
opposite spins and the Stoner parameter / is given as 

/ = UN. (8) 

Eq.Elmakes it clear that, for a finite /, singularities in 
the DOS can be expected to induce spontaneous magne- 
tization 0. In QDs with high symmetry (e.g., spherical 
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FIG. 2: Typical EMD in a QD, n(p), and the magnitude of 
its first derivative, |n'(p)[. The position of the peak in |n'(p)| 
defines the QD "Fermi momentum", while its full-width-at- 
half-maximum defines Ap. 



FIG. 3: Ap in reduced units of pf (solid line), and spin 
polarization £ (dashed) vs. QD radius 7?. Simulations are 
based on the model Hamiltonian of Eq. Q using parameters 
discussed in the text. 



or cubic dots), symmetry related degeneracies will gener- 
ally enhance DOS peaks. On the other hand, symmetry 
breaking effects and disorder will smear out DOS 
peaks and reduce the tendency for the system to magne- 
tize spontaneously. 

We discuss next the EMD with an eye towards identi- 
fying signatures of spin polarization in a QD. The EMD 
is defined by 

/oc 
dLjf(u)A(p,u) (9) 
-oo 

where / is the Fermi function and 

A(p, u) = -2ImG fl (p, to), (10) 

is the spectral function. The one particle Green's func- 
tion G R (p,ui) and its imaginary part can be evaluated 
exactly for the model Hamiltonian of Eq. ^ The typical 
behavior of the EMD and its derivative is shown in Fig. [2] 
The region of rapid variation in n(p) can be characterized 
via the position, pp, of the peak in |n'(p)| and the associ- 
ated full-width-at-half-maximum, Ap. In the bulk limit 
in a metallic system, the EMD in general contains Fermi 
surface (FS) breaks in the first Brillouin zone (BZ) and 
at the Umklapp images of the FS in higher BZs. Cor- 
respondingly, the first derivative of the EMD develops 
(5-function peaks. Although in a finite system there can- 
not be breaks in the EMD, we may nevertheless refer to 
Pf loosely as the QD "Fermi momentum" for simplicity. 

Fig. shows the simulated behavior of Ap in a QD as 
a function of the radius R at a fixed electron density. Ap 
is seen to display peaks, which are well correlated with 



those in the magnetization £. The reason for this corre- 
lation between Ap and C is that in the polarized system, 
in effect, there are two separate momentum distributions 
for the up and down spin electrons. The two associated 
" Fermi" momenta then give rise to two different peaks in 
\n'(jj)\, which appear as increased broadening of Ap in 
the total momentum density of the interacting system. 
In fact the R dependence of Ap can be fitted as [4J 

Ap/p F = Ap /p F + cA(C-Co)i?, (11) 

where Co is the spin polarization in the weakly interact- 
ing case (U — > 0), and c is a constant, which depends on 
various QD parameters. Eq. 1111 can be used to extract 
the polarization £ from the measured R-dependence of 
Ap. These considerations indicate that peaks in Ap pro- 
vide a distinct signature of polarization of a QD and that 
spectroscopies sensitive to momentum density can play a 
useful role in this connection [1J|. These results argue 
for investigations of the QDs using Compton scattering 
and positron annihilation experiments. 

It should perhaps be noted that Compton scattering 
and positron annihilation spectroscopies have developed 
by now into standard probes of the EMD in materials. 
Recent positron annihilation measurements on CdSe QDs 
show that the effect of quantum confinement results in 
an increased width Ap of the Fermi edge in momentum 
space. The width Ap seems to follow an inverse square 
law 1 /d 2 with particle diameter d 0, ^| , in contrast 
to the 1/d law expected for the confined homogeneous 
electron gas |l7j . 



4 



III. A QUANTUM MONTE CARLO 
CALCULATION 



The discussion of the preceding section is based on the 
Hamiltonian of Eq. ^and it is thus limited by the form of 
the Hartree-Fock many-body wavefunction given by Eq. 
|5J We now consider a more general many-body wave- 
function and apply the VQMC approach to focus 
particularly on understanding the nature of the effective 
interaction between electrons of the same and opposite 
spins in a QD. The specific wavefunction used is 



0.7 



where 



* = JD^D 1 , 



J = cxp(2_^ Uij) 

Kj 



(12) 



(13) 



is the Jastrow factor, which is expressed in terms of a 
product involving two-body correlation factors, tty. The 
must fulfill the so-called cusp condition, i.e. the sin- 
gularities associated with the kinetic energy must cancel 
those arising from the Coulomb potential in the micro- 
scopic Hamiltonian. A simple form is given by (lif 



u(r) 



l + /3r 



(14) 



where f3 is a variational parameter. When J = 1 or 
[3 — > oo, ^ reduces to the form of an unrestricted Hartree- 
Fock wave function. 

The average repulsion energy between two electrons 
can be expressed as an integral of the spin resolved radial 
pair correlation function 



9*,*' ( r ) = ^XX^<Ws / s ( 



ry)|*(0rdR 
(15) 

where R = (ri, ...,rjv), £ = (R,eri, ...,<tn), r« = l^-r^l, 
5 a ,cr' is the spin projector, and n is the average electron 
density. The pair interaction energy between two elec- 
trons of like spins is then given by 



1 



nga,<r (r) 



(No - 1) 

and for electrons of opposite spins is 

1 f ng„- a (r) 



dr. 



Ve,-a = 



N„ 



dr. 



(16) 



(17) 



The energy penalty U for having a pair of electrons with 
opposite spins can therefore be obtained from the average 



(18) 



As an example, we have investigated a QD consisting 
of 12 electrons enclosed in a 2D square well of size I = n 



Opposite spin 

Same spin 




FIG. 4: Radial pair correlation functions ga,a{r) (dashed line) 
and g CTj _cr(r) (solid line) for 12 electrons enclosed in a 2D 
square well of size I = it a* B . 



a* B [19j. Here and elsewhere in this section, it is con- 
venient to use the modified atomic units a* B for length 
and H* for energy, which are renormalized atomic units 
obtained from the effective electron band mass m* and 
the dielectric constant of the material e [2(J ■ 

We describe the confinement in the xy plane by an 
infinite hard-wall potential, therefore the single-electron 
states in the square QD are 



, a (x,y) = -s'm(n x x) sm(n y y). 

7T 



(19) 



For N = 12, one obtains a closed shell system (iVf = JVj.) 
with zero net magnetization. 

In the Hartree-Fock limit ((3 — ► oo), a Monte Carlo cal- 
culation gives the total energy of the 12-electron QD to be 
107.785 ±0.002 H*. The parameter (3 was then optimized 
via the Stochastic Gradient Approximation (SGA) (lij| . 
In the SGA, at each step n, the value of an observable x 
is updated with a recursive calculation of the mean: 



Xn—l 



-(X n -1 

n 



Xji) 



(20) 



At the optimal f3 = 1.53, the total energy is found to 
be 103.237 ±0.001 H* with \dE/d(3\ < 10~ 3 Wa* B . Us- 
ing this value of (3, the results of a relatively noiseless 
calculation of the spin-dependent pair correlation func- 
tions g a <7 r are shown in Fig. 21 As expected, the dis- 
tribution ga t<7 (r) (dashed line) is seen to vanish at r = 0, 
reflecting the presence of an " exchange hole" surrounding 
like spins due to the Pauli exclusion principle. Electrons 
of unlike spins, on the other hand, still tend to avoid 
each other due to Coulomb repulsion, which induces a 
"Coulomb" or "correlation" hole in g a ^ a (r) (solid line). 
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Fig. 0] shows that the "hole," or the region of depleted 
electron density, excludes a larger number of electrons 
and extends to a larger distance for like spins than for 
unlike spins. The decrease of g a . a > after r rs 1 a* B is a 
geometrical effect due to the finite size of the QD. The 
use of these correlation functions in Eqs. 16-18 yields, 
V n = 1.039 H*, V\i = 0.760 H*, and U = 0.279 H*. 
The corresponding values in the Hartree-Fock limit are: 
V n = 1.176 H*, V n = 0.848 H*, and U HF = 0.329 
H*. The Jastrow wavefunction thus leads to a reduc- 
tion in the energy penalty for creating a pair of opposite 
spins. The overall effect however is relatively small in 
that the effective U for the Jastrow wavefunction is only 
15% smaller than Uhf- These results suggest that, de- 
spite its simplicity, the model Hamiltonian of Eq. ^ is 
capable of providing a reasonable description of the elec- 
tron gas in QDs. Moreover, we have explicitly verified 
the oscillations of the magnetization in the 3D-spherical 
quantum dots containing up to 8 electrons by perform- 
ing VQMC simulations with the SGA optimization of 
the total energy. The oscillations in magnetization with 
QD radius predicted on the basis of this Hamiltonian are 
presumably robust to electron correlation effects missing 
implicitly in the Hartree-Fock form of its solution. 



IV. SUMMARY AND CONCLUSIONS 

We discuss issues of spin polarization and momentum 
density with focus on the treatment of correlation ef- 
fects in electron gas confined within a quantum dot. In 
this connection, we first consider selected results based 
on a relatively simple model Hamiltonian of Ref. |4| in 
which interactions are restricted via a parameter U to 
be non-zero only for electrons of opposite spins. This 
model Hamiltonian is solvable exactly and admits a so- 
lution of the Hartree-Fock form. Moreover, it displays re- 
markable oscillations in spin polarization with dot radius, 
which leave distinct signature in the electron momentum 
density. In order to gain insight into correlation effects 
more generally, we have carried out VQMC calculations 
on a square dot containing 12 electrons using a Jastrow- 
Slater form of the many body wavefunction. The effec- 
tive U value for the Jastrow-Slater wavefunction is found 
to be only 15% smaller than for the Hartree-Fock case. 
On the whole, we conclude that spectroscopies sensitive 
to electron momentum density-Compton scattering and 
positron annihilation in particular-can potentially help 
delineate spin polarization effects in quantum dots. 

This work is supported by the US Department of En- 
ergy contract DE-AC03-76SF00098, and benefited from 
the allocation of supercomputer time at the NERSC and 
the Northeastern University's Advanced Scientific Com- 
putation Center (ASCC). 



[1] L. P. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. 

Tarucha, R. M. Westervelt, and N. S. Wingreen, in Meso- 

scopic Electron Transport, NATO Science Series: E, Vol. 

345, edited by L. L. Sohn, L. P. Kouwenhoven, G. Schon, 

(Kluwer, Dordrecht, 1997). 
[2] S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 

1283 (2002). 

[3] H. Jiang, H.U. Baranger, and W. Yang, Phys. Rev. Lett. 

90, 026806 (2003). 
[4] R. Saniz, B. Barbiellini, A. B. Denison, A. Bansil, Phys. 

Rev. B 68, 165326 (2003). 
[5] I. I. Yakimenko, A. M. Bychkov, and K.-F. Berggren, 

Phys. Rev. B 63, 165309 (2001). 
[6] N. Zabala, M. J. Puska, and R. M. Nieminen, Phys. Rev. 

Lett. 80, 3336 (1998) and Phys. Rev. B 59, 12652 (1999). 
[7] G. L. Gutsev, S. N. Khanna, and P. Jena, Phys. Rev. B 

62, 1604 (2000); for a general reference see also Quantum 

Phenomena in Clusters and Nano structures, Ed. S. N. 

Khanna and A. W. Castleman, Jr., (Springer, New York, 

2003). 

[8] D. J. Reilly, T. M. Buehler, J. L. OBrien, A. R. Hamilton, 
A. S. Dzurak, R. G. Clark, B. E. Kane, L. N. Pfeiffer, and 
K. W. West Phys. Rev. Lett. 89, 246801. 

[9] W. H. Li, S. Y. Wu, C. C. Yang, S. K. Lai, K. C. Lee, H. 
L. Huang and H. D. Yang, Phys. Rev. Lett. 89, 135504 

(2002) ; H.Hori, T. Teranishi, Y. Nakae, Y. Seino, M. 
Miyake and S. Yamada, Phys. Lett. A 263, 406 (1999); H. 
Hori, Y. Yamamoto, T. Iwamoto, T. Miura, T. Teranishi, 
M. Miyake, cond-mat/0306261 to appear in Phys. Rev. 
B. 

[10] The Coulomb matix elements and therefore U in general 

depend on the size of the system. 
[11] O. Gunnarson, J. Phys. F: Metal Phys. 6, 587 (1976). 
[12] CM. Hurd, Electrons in Metals, Krieger (Malabar, 

Forida, 1981). 

[13] Ph. Jacquod and A. D. Stone, Phys. Rev. Lett. 84, 3938 
(2000) 

[14] Complex magnetic structures may be expected in dots of 
non-symmetric shapes, including the possibility of spiral 
spin arrangments and the formation of domains for dots 
of large sizes. It will be interesting to consider associated 
signatures in the momentum density. Magnetic Compton 
scattering with its ability to focus on specific magnetic 
states may prove to be a particularly useful spectroscopy 
in this regard (see, e.g., P. A. Montano, Yinwan Li, J. 
F. Mitchell, B. Barbiellini, P.E. Mijnarends, S. Kaprzyk, 
and A. Bansil, J. Phys. Chem. Solids, this volume). 

[15] M.H. Weber, K.G Lynn, B. Barbiellini, PA. Sterne, 
A.B. Denison, Phys. Rev. B 66R 041305 (2002). 

[16] A.B. Denison, R. Meulenberg, S.W.H. Eijt, A. Van Veen, 
P.E. Mijnarends, B. Barbiellini, A. Bansil, C. Fischer, 
M.H. Weber, K.G. Lynn, Mat. Sci. Forum 445-446, 395 
(2004). 

[17] R. Saniz, B. Barbiellini and A. Denison, Phys. Rev. B 

65, 245310 (2002). 
[18] A. Harju, B. Barbiellini, S. Siljamaki, R.M. Nieminen 

and G. Ortiz, Phys. Rev. Lett. 79, 1173 (1997). 
[19] E. Rasanen, H. Saarikoski, V.N. Stavrou, A. Harju, M.J. 

Puska and R.M. Nieminen, Phys. Rev. B 67, 235307 

(2003) . 

[20] The band mass m* and the dielectric constant e of the 



6 



polarizable background influence the modified atomic stance, in GaAs, using m* = 0.067m and e = 12.4, yields 

units: H*= (e 2 /m*) H, a B — (m*/e)as, where H is the H*= 11.8572 meV and a* B = 9.79 nm. 

Hartree energy unit and as is the Bohr radius. For in- 



